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Gaussian random-matrix ensembles defined over the tangent spaces of the large families of Car- 
tan’s symmetric spaces are considered. Such ensembles play a central role in mesoscopic physics 
as they describe the universal ergodic limit of disordered and chaotic single-particle systems. The 
generating function for the spectral correlations of each ensemble is reduced to an integral over a 
Riemannian symmetric superspace in the limit of large matrix dimension. Such a space is defined 
as a pair (G/H,M I ) where G/H is a complex-analytic graded manifold homogeneous with respect 
to the action of a complex Lie supergroup G, and M r is a maximal Riemannian submanifold of the 
support of G/H. 
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I. INTRODUCTION 

The mathematics of supersymmetry, though conceived and developed in elementary particle theory, has been applied 
extensively to the physics of disordered metals during the past decade. Improving on earlier work by Wegner 
Efetov j3j] showed how to approximately map the problem of calculating disorder averages of products of the energy 
Green’s functions for a single electron in a random potential, on a supersymmetric nonlinear a model. Later it was 
shown 0 that the same nonlinear a model describes the large-TV limit of a random-matrix ensemble of the Wigner- 
Dyson [g type. Since then, Efetov’s method has evolved into a prime analytical tool in the theory of disordered or 
chaotic mesoscopic single-particle systems. Competing methods are limited either to the diffusive regime (the impurity 
diagram technique), or to isolated systems in the ergodic regime (the Dyson-Mehta orthogonal polynomial method), 
or to quasi-one-dimensional systems (the DMPK equation). In contrast, Efetov’s method is applicable to isolated 
and to open systems in the diffusive, ergodic, localized, and even ballistic regime, to both spectral correlations and 
transport properties, and it can in principle be used in any dimension. This versatility has engendered a large body 
of nontrivial applications, many of which are outside the range of other methods. Of these, let me mention: (i) the 
Anderson transition on a Bethe lattice [[i] |8| , (ii) localization in disordered wires jii] [l3| , (iii) multifractality of energy 
eigenstates in two dimensions jbi] [n|, (iv) weak localization and conductance fluctuations of chaotic billiards strongly 
coupled to a small number of scattering channels [jl7],|l8j and, most recently, (v) a theoretical physicist’s proof of the 
Bohigas-Giannoni-Schmit conjecture for chaotic Hamiltonian systems [ pT)||20| l . 

In spite of these manifest successes, Efetov’s supersymmetry method has been ignored (for all that I know) by 
mathematical physicists. This is rather unfortunate for several reasons. First, an infusion of mathematical expertise 
is needed to sort out some matters of principle and promote the method to a rigorous tool. Second, various extensions 
of currently available results seem possible but have been hindered by the lack of mathematical training on the part 
of the condensed matter theorists applying the method. And third, the geometric structures underlying Efetov’s 
nonlinear a models are of exquisite beauty and deserve to be studied in their own right. Part of the reason why 
neither mathematicians nor mathematical physicists have monitored or contributed to the development, may be that 
there does not exist a concise status report that would appeal to a mind striving for clarity and rigor. Hence the first, 
and very ambitious, motivation for getting started on the present paper was to make an attempt and partially fill the 
gap. 

Another objective is to report on a recent extension of the supersymmetry method to random-matrix theories 
beyond the standard Wigner-Dyson ones. In her study of Anderson localization in the presence of an A-B sublattice 
symmetry, Gade pl| noticed that the manifold of the nonlinear cr model is promoted to a larger manifold at zero 
energy. The same phenomenon occurs in the chiral limit of the QCD Dirac operator at zero virtuality |^2| . For 
several years it remained unclear how to handle this enlargement of the manifold in the supersymmetric scheme. 
(Gade used the replica trick instead of supersymmetry.) The key to solving the problem can be found in a paper 
by Andreev, Simons, and Taniguchi |2f| who observed that what one needs to do is to avoid complex conjugation 
of the anticommuting variables. In the present paper I will elaborate on this observation and cast it in a concise 
mathematical language. Moreover, I will show that the same technical innovation allows to treat the random-matrix 
theories that arose j24],|25|] in the stochastic modeling of mesoscopic metallic systems in contact with a superconductor. 

An outline of the basic mathematical structure is as follows. Consider a homogeneous space G/H, where G and H 
are complex Lie supergroups, and regard G/H as a complex-analytic (p, gj-dimensional supermanifold in the sense of 
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Berezin-Kostant-Leites To integrate its holomorphic sections, select a closed, oriented and realp-manifold M r 

contained in the support M = Go/Hq of the supermanifold. The natural (invariant) supergeometry of G/H induces 
a geometry on M r by restriction. If this geometry is Riemann and M r is a symmetric space, the pair ( G/H,M r ) is 
called a Riemannian symmetric superspace. This definition will be shown to be the one needed for the extension 
of the supersymmetry method beyond Wigner-Dyson. The difficulties disordered single-particle theorists had been 
battling with were caused by the fact that the exact sequence 

0 —> nilpotents —> G/H —> M —> 0 


does not, in general, reduce to an exact sequence of sheaves of real-analytic sections terminating at the Riemannian 
submanifold M r . 

When integrating the invariant holomorphic Berezin s uper form on G/H , one must pay careful attention to its 
coordinate ambiguity. This subtle point is reviewed in Sec. II A. After a brief reminder of the procedure of Grassmann- 
analytic continuation (in Sec. (IB), the complex Lie supergroups Gl(m|n) and Osp(m|2n) (in Sec. [IIC| ), and Car tan’s 
symmetric spaces (in Sec. EH. the details of the definition of Riemannian symmetric superspaces are given in Sec. |ll F| . 
Table 2 lists the large families of these spaces. 

Sec. HI, the largest of the paper, treats the Gaussian random-matrix ense mble defined over the symplectic Lie 
algebra sp(IV), by an adaptation of Efetov’s method. A simple example (Sec. Ill A) illustrates the general str ategy. 
Details of the method, including a complete justification of all manipulations involved, are presented in Secs. Ill B 

|III F . Theorem 3.3 expresses the Gaussian ensemble average of a product of n ratios of spectral determinants 
as a superintegral. Theorem 3.4 reduces this expression to an integral over the Riemannian symmetric superspace 
Osp(2n|2n)/Gl(n|n) with M r = (SO*(2n)/U(n)) x (Sp(n)/U(n)), in the limit N —> oo. 

According to Cartan’s list, there exists eleven large families of symmetric spaces. Ten of these correspond to 
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universality classes that are known to describe disordered single-particle systems in the ergodic regime 124 2E 
class singled out for detailed treatment in Sec. Ill describes mesoscopic normal-superconducting hybrid systems 
with time-reversal symmetry broken by a weak magnetic field. The remaining nine classes are briefly discussed in 
Sec. IV. Each of them is related, by the supersymmetry method, to one of the large families of Riemannian symmetric 
superspaces of Table 2. A summary is given in Sec. |V|. 


II. RIEMANNIAN SYMMETRIC SUPERSPACES 


A. The Berezin integral on analytic supermanifolds 


Let A(U) denote the algebra of analytic functions on an open subset U of p-dimensional real space. By taking 
the tensor product with the Grassmann algebra with q generators one obtains A(U) (g> A(R 9 ), the algebra of analytic 
functions on U with values in A(R 9 ). Multiplication on A(R 9 ) is the exterior one, so the algebra is supercommutative 
(or graded commutative). The object at hand serves as a model for what is called a real-analytic (p, gj-dimensional 
supermanifold (or graded manifold |28|) in the sense of Berezin, Kostant, and Leites (BKL) | [fti[|27| ] : which, precisely 
speaking, is a sheaf of supercommutative algebras A with an ideal J\f (the nilpotents), such that M ~ A/M is 
an analytic p-manifold and on a domain U C M, A splits as A(U) (g> A(R 9 ). The global sections of the bundle 
A —> M are called superfunctions, or functions for short. M is called the underlying space, or base, or support, of 
the supermanifold. M will be assumed to be orientable and closed (dM = 0). 

The calculus on analytic supermanifolds is a natural extension of the calculus on analytic manifolds. Functions are 
locally expressed in terms of (super-)coordinates (x;£) := (x 1 ,... ,x p ;£ x , ...,£ 9 ) where x l (£ J ) are even (resp. odd) local 
sections of A. If (x;£) and ( 3 / 577 ) are two sets of local coordinates on domains that overlap, the transition functions 
y 1 = f l (x; £) and p - 3 = qA(x; £) are analytic functions of their arguments and are consistent with the Z 2 -grading of A. 

In what follows the focus is on the theory of integration on analytic supermanifolds. Recall that on p-manifolds the 
objects one integrates are p-forms and their transformation law is given by 

d y 1 A ... A d y p = dx 1 A ... A dx p Det (dy l /d, ad) . 

The obvious (super-)generalization of the Jacobian Det {dy l /dx^) is the Berezinian Q 


Ber (y, 77 /x, £) := SDet 


f dy l /dxi 

ydrii/dx 0 


dyt/dg 

drji/d\ j 


where SDet is the symbol for superdeterminant. Guided by analogy, one postulates that an integral superform ought 
to be an object D transforming according to the law 
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D(y, rf) = D(x , £) Ber (y, rj/x, £) . 


Cl) 


A natural candidate would seem to be 


D{x,£) := drr 1 A ... A dx p g> 

which is a linear differential operator taking superfunctions / into p -forms D[f ] (<9ji denotes the partial derivative 
with respect to the anticommuting coordinate £ l ). The p-form D[f} can be integrated in the usual sense to produce 
a number. However, the transformation law for D( a;,£) turns out to be not quite (0) but rather 

D(y,v) = D(x,£) Ber (y, r?/x, 0 + 0. ( 2 ) 

An explicit description of the term /? on the right-hand side, here referred to as the anomaly , was first given by 
Rothstein |iot . It is nonzero whenever some even coordinate functions are shifted by nilpotent terms. Its main 
characteristic is that on applying it to a superfunction /, one gets a p-form that is exact: /3[f] = d(a[/]). 

The existence of an anomaly in the transformation law for D{x,£) leads one to consider a larger class of objects, 
namely A P (M) <8> .4 2?, the sheaf of linear differential operators on A with values in the p-forms on M. (A P (M) <g> jJ) 
naturally is a right A-module.) To rescue the simple transformation law (|l[) one usually passes from A P (M) (g> XD to 
its quotient by the anomalies po|| . In order for the integral to be well-defined over the quotient, one must take the 
functions one integrates to be compactly supported. 

Sadly, this last option is not available to us. The functions that will be encountered in the applications worked out 
below, do not ever have compact support but are analytic functions instead. When integrating such functions, we 
need to work with the full transformation law (^|), which includes the anomaly. 

Another way of avoiding the anomaly is to arrange for the transition functions never to shift the even coordinates 
by nilpotents, by constructing a restricted subatlas |3l]]. However, because the concept of a restricted subatlas is 
somewhat contrived, this approach has been found to be of limited use in the type of problem that is of interest here. 

To arrive at a definition of superintegration that is useful in practice, we proceed as follows. The supermanifold is 
covered by a set of charts with domains [/; and coordinates (x/j \, (i = 1,..., n ). On chart i let u)i := D(x 

with u>i a local section of A, and let a* £ A p- 1 (M) ® XD\u t - Partition M into a number of consistently oriented p-cells 
D i,..., D n , with Di contained in C/j. For i < j put Py := dDi D dDj and, if D; Lj is nonempty and is a (p — l)-cell, fix 
its orientation by dDi = +Py + ... . 

Definition 2.1: A collection is called a Berezin measure to if the conditions 

Ui/ujj = Ber (i/j), (3) 

cji + da j = u)j + d aj , (4) 

are satisfied on overlapping domains. The Berezin integral / i—> f M lu [/] is defined 


>\f\ = 


! M 


E 

2=1 ' 


A[/] 


Di 


i<j ' 


iAf] 


Di 


( 5 ) 


where a^- = a^ — aj. The quantities u>i and a,; are called the principal term and the anomaly of the Berezin measure 
on chart i. 

Remark 2.2: The conditions (|j) and ([|) ensure the existence of a global section to £ A P (M) (g) whose local 

expression in chart i is uji + da^. The existence of to means that the distribution (|5|) is independent of the coordinate 
systems and the cell partition chosen. Because (^) depends only on the differences a, — aj, one can gauge the anomaly 
to zero on one of the charts without changing the Berezin integral. 

Example 2.3: Consider the real supersphere S p l 2 , a (p, 2)-dimensional supermanifold with support S p , which is the 
space of solutions in (p + 1, 2) dimensions of the quadratic equation 

*o+^i+ ••• + Xp + 2|i6 = 1- 

Cover S p by two domains 1 and 2 obtained by removing the south (xq = —1) or north pole (xq = +1). Introduce 
stereographic coordinates (x\,..., x p \ £i, £ 2 ) and (y±,..., y p \ iji, 772 ) for S p l 2 on these domains with transition functions 

pi = —xi/R 2 , yi = Xi/R 2 (i = 2, ...,p), rjj = £j/R 2 (j = 1,2) 
where R 2 = x i + 2^i^ 2- (The minus sign preserves the orientation.) Consider 
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Wi =D(x,t)o (l + J2 x i + 2 66) P+2 , 
w 2 = D( y , rj) o (1 + Yjyf + P+2 , 

fV r 2 'l (p_2) / 2 

“ 12 = ~ n (l+E^)H 0 ° ^ 2 ’ 

where fl = (^ . a; 2 ) _p / 2 ^? =1 (—l) l dxi A ... Ad:Cj_i /\Xi&Xi + i A ... Adx p is the solid-angle (p— l)-form in p dimensions. 
It is not difficult to check by direct calculation that u>i, u >2 and au = aj. 2 , c *2 = 0 obey the relations © and d). 
Hence, they express a glo bally defined Berezin measure u> in the sense of Definition 2.1. (The geometric meaning of 
u> will be specified in Sec. |TQ.) For p > 3, the anomaly a \2 scales to zero when ^xf —> oo, so we may shrink cell 2 
to a single point (a set of measure zero) and compute the Berezin integral simply from 

[ v[f]= [ D(x,£) (l + J2x 2 + 2£i&y P+2 f(x;£). 

J SP Jrp 

In these cases we can get away with using only a single chart. The situation is different for p = 2 and p = 1. In the 

first case the anomaly is scale-invariant (the solid angle is) and by again shrinking cell 2 to one point (the south pole 

(*,») = (0,0) on /) we get 



D{x,£)f(x-,£) +4tt/ 


south pole 


In particular, J g2 w[l] = 47r. For p = 1 the anomaly diverges at x = 0 and x = oo. In this case the general formula 
(j5|) must be used, and one finds / gl w[l] = 0. 


B. Grassmann-analytic continuation 


In the formulation of BKL, the vector fields of a supermanifold do not constitute a module over A but are constrained 
to be even derivations of A. which is to say that their coordinate expression is of the form 


x = r(x;()^ + V ‘M^ 


where f l and are even and odd superfunctions respectively. Unfortunately, this formulation is too narrow for 
most purposes. The reason is that in applications one typically deals not with a single supermanifold but with 
many copies thereof (one per lattice site of a lattice-regularized field theory, for example). So in addition to the 
anticommuting coordinates of the one supermanifold that is singled out for special consideration, there exist many 
more anticommuting variables associated with the other copies of the supermanifold. When the focus is on one 
supermanifold, these can be considered as “parameters”. Often one wants to make parameter-dependent coordinate 
transformations, leading to coefficients /*!...i n (a;) in the expansion /(a;;£) = fh~-i n (a;)£ n ■■■C n that depend on 
extraneous Grassmann parameters. (For example, when the supermanifold is a Lie supergroup, it is natural to 
consider making left and right translations g i—> gLggR-) The upshot is that one wants to take A as a sheaf of graded 
commutative algebras not over R but over some (large) parameter Grassmann algebra A (the Grassmann algebra 
generated by the anticommuting coordinates of the “other” supermanifolds). Making this extension, which is called 
“Grassmann-analytic continuation” in J|i| , one is led to consider the more general class of vector fields of the form 

X = f i (x,Z-,f3)-^+y(x,t;-,P)-^ ( 6 ) 

where the symbol 3 stands for the extra Grassmann parameters and the dependences on these are such that f‘ and 
continue to be even and odd respectively (the Z 2 -grading of A after Grassmann-analytic continuation is the natural 
one). 

The vector fields (;(i ) still are even derivations of the extended algebra. One can go further by demanding that Der.4 
be free over A and including the odd ones, too. When that development is followed to its logical conclusion, one 
arrives at Rothstein’s axiomatic definition |j^] of supermanifolds, superseding an earlier attempt by Rogers ]33| , [34[ . 
Although there is no denying the elegance and consistency of Rothstein’s formulation, we are not going to embrace it 
here, the main reason being that odd derivations will not really be needed. For the purposes of the present paper we 
will get away with considering vector fields of the constrained form (^|). 
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C. The complex Lie supergroups Gl(m|n) and Osp(m|2n) 

The supermanifolds we will encounter all derive from the complex Lie supergroups |^9|,^8| Gl(m|n) and Osp(77i|2n), 
by forming cosets. The definition of Gl(m|?z) rests on the notion of an invertible supermatrix g = ( 900 901 ) where 

V ffio Sii/ 

goo , goi, <710 and <7n are matrices of size m x m, m x n, n x to, and n x n. The supermanifold structure of Gl(m|n) 
comes from taking the matrix elements of goo and gn (goi and gio) for the even (resp. odd) coordinates on suitable 
domains of the base M = Gl(m, C) x Gl(n, C). The Lie supergroup structure derives from the usual law of matrix 
multiplication. 

For m / n, it is common practice to split, off from Gl(m|n) the Gl(l)-ideal generated by the unit matrix, so as 
to have an irreducible Lie superalgebra [|35[ , p6| . For m = n, which turns out to be the case of most interest here, 
one ends up having to remove two Gl(l)’s, one generated by the unit matrix and the other one by the superparity 
matrix a = diag(l„, — l n ). And even then the Lie superalgebra is not irreducible in a sense, for the Killing form 
STr ad(X)ad(K) vanishes identically. We therefore prefer to take Gl(m|n) as it stands (with no ideals removed) 
and replace the Killing form by the invariant quadratic form B(X, Y) = STrXK, which is nondegenerate in all cases 
(including m = n). 

The complex orthosymplectic Lie supergroup Osp(m|2n) is defined as a connected subgroup of Gl(m|2n) fixed by 
an involutory automorphism q <—> f(q) = to -1 t -1 , where r is supersymmetric (r = r T cr = crT T )Jd The support of 
Osp(ro|2n) is SO(m,C) x Sp(ra,C). ' 

The action of a Lie supergroup on itself by left and right translations gives rise to right- and left-invariant vector 
fields. A Berezin measure on a Lie supergroup is said to be invariant, and is called a Berezin-Haar measure, if its Lie 
derivatives |2^] with respect to the invariant vector fields vanish. 

Given a Lie supergroup G and a subgroup H , the coset superspace G/H is defined by decreeing that the structure 
sheaf of the coset superspace is a quotient of sheaves. The action of G on G/H by left translation gives rise to so-called 
Killing vector fields. A Berezin measure on G/H is called invariant if its Lie derivatives with respect to the Killing 
vector fields are zero. 

If Osp R (?n|2?i) denotes the orthosymplectic supergroup over the reals, the coset space Osp r (to+ l|2?r)/Osp R (m|2n) 
can be identified with the real supersphere S m l 2ra . The Berezin measure discussed in Example 2.3 is invariant with 
respect to the action of Osp R (p+ 112) on S^ 2 and can be viewed as the “volume superform” of S p l 2 . Hence we can 
restate the results of that example as follows: vol(S 2 l 2 ) := / g2 w[l] = 47 t and vo^S 1 * 2 ) := / gl w[l] = 0. 


D. Holomorphic Berezin measures on complex-analytic supermanifolds 


To go from real-analytic supermanifolds to complex-analytic ones, one replaces the structure sheaf A by a sheaf 
of graded commutative algebras H over C such that M ~ Ti/Af is a complex manifold and H. is locally modeled by 
H(U) <g> A(C 9 ), where H(U) is the algebra of holomorphic functions on U C M. The natural objects to consider then 
are holomorphic superfunctions, i.e. global sections of the bundle H —► M. In local coordinates z 1 , ...,z p ]C, 1 , ...,C 9 
such sections are written f(z', £). Grassmann-analytic continuation is done as before when needed. A Berezin measure 
on a complex-analytic ( p , q)-dimensional supermanifold is a linear differential operator w th at takes holomorphic 
superfunctions / into holomorphic p -forms u>[f] on M. The statements made in Sec. IIA about the anomalous 
transformation behavior of Berezin measures apply here, too (mutatis mutandi). 

To define Berezin’s integral in the present context, one more piece of data must be supplied, namely a real 75 - 
dimensional submanifold M r C M over which the holomorphic 75 -form ui[f] can be integrated to produce a complex 
number. Thus, given to and M r , Berezin’s integral is the distribution 


/' 


I M t 


<*>[/]■ 


(7) 


Let me digress and mention that this definition, natural and simple as it is, was not “discovered” by the random- 
matrix and mesoscopic physics community (including myself) until quite recently. With one notable exception ||23|j , all 
past superanalytic work on disordered single-particle systems employed some operation of “complex conjugation” of 
the Grassmann generators - namely an adjoint of the first or second kind [p9| - to make the treatment of the ordinary 


1 Supertransposition T is an operation with the properties (AB) T = B T A T and A TT = aAa. 


5 




(“bosonic”) and anticommuting (“fermionic”) degrees of freedom look as much alike as possible. Presumably this 
was done because it was felt that such egalitarian treatment is what is required by the principle of “supersymmetry”. 
Specifically, a reality constraint was imposed not just on the underlying space M (fixing M r ) but on the entire 
structure sheaf to reduce Ti to a sheaf of algebras over R. Although this reduction can be done with impunity in some 
cases (namely the classic Wigner-Dyson symmetry classes), it has turned out to lead to insurmountable difficulties 
in others (the chiral and normal-superconducting symmetry classes). A major incentive of the present paper is to 
demonstrate that the construction (j?j) is in fact the “good” one to use for the application of supermanifold theory to 
disordered single-particle systems in general. Although that construction may hurt the physicists’ aesthetic sense by 
“torturing supersymmetry”, it should be clear that we are not breaking any rules. Recall that according to Berezin, 
superintegration is a two-step process: first, the Fermi-integral (i.e. differentiation with respect to the anticommuting 
coordinates) is carried out, and it is only afterwards that the ordinary (Bose-) integrals are done. When the sequential 
nature of the Berezin integral is taken seriously, there is no compelling reason why one should ever want to “complex 
conjugate” a Grassmann variable. In the present paper, we take the radical step of abandoning complex conjugation 
of Grassmann variables altogether. 

Example 2.f: The simplest nontrivial example [^3| is given by Gl(l|l), the Lie supergroup of regular complex 
2x2 supermatrices g = ^ ^ with support M = G1(1,C) x G1(1,C). The Berezin-Haar measure on Gl(l|l) is 

u> = (2ni)~ 1 D(ad; /Ty) where D(ad;/3 7 ) = da A dd ® dpd 7 . Solving the regularity conditions a ^ 0 and d / 0 by 
parameterizing Gl(lll) through its Lie algebra, g = exp ( 21 1 ), one finds 

2 "“ = D(Z ' Z2:<1<2 > ° (1 - e - ( -»Xe«- - 1) 

- (-»"(■“ - ■-> - ,(!■;, • «, m 

Note that this expression is holomorphic in a neighborhood of the origin z± = z 2 = 0. The first term on the right-hand 
side is the principal term, and the second one is the anomaly of to in these coordinates. To integrate to, one might be 
tempted to choose for M r the U(l) x U(l) subgroup defined by Re(zi) = 0 = Re(z 2 ). However, since the rank-two 
tensor STr dgdg -1 = dada -1 — dddc?” 1 + nilpotents = —dzf + d z\ + ... is not Riemann on U(l) x U(l), this will not 
be the best choice. A Riemannian structure is obtained by taking M r = M + x S 1 defined by Im( 2 i) = 0 = Re(z 2 ). To 
compute / R + xg i ai[f] we may use a single cell 

D : —00 < x < + 00 , — 7 r < y < +7T, 

where x = Re(zi) and y = Im(^ 2 ). The boundary 9D consists of the two lines y = —n and y = n (x € R). Using (||), 
paying attention to the orientation of the boundary, and simplifying terms, one finds the following explicit expression 
for the integral of ui: 



By construction, this Berezin integral is invariant under left and right translations f(g) 1 —> f{gLggn)- Evaluation 
gives / R+xg 1 uj[ 1] =1^0. The naive guess would have been J w[l] = (2ni)~ 1 f da Add dpd 7 = 0 due to dgd-y -1 = 0. 
Such reasoning is false because f R+ da = 00 . 

E. Symmetric spaces: a reminder 

A Riemannian (globally) symmetric space is a Riemannian manifold M such that every p £ M is an isolated fixed 
point of an involutive isometry. (In normal coordinates x l centered around p , this isometry is given by x l 1 —> — x l .) 
This definition implies (cf. |p7| ) that the Riemann curvature tensor is covariantly constant, so that “the geometry 
is the same everywhere”. The curvature can be positive, negative or zero, and the symmetric space is said to be of 
compact, noncompact or Euclidean type correspondingly. 
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According to Cartan’s complete classification scheme, there exist tenQ large classes of symmetric spaces. Apart 
from some minor modifications the motivation for which is given presently, these are the entries of Table 1: 


class 

noncompact type 

compact type 

A 

AI 

All 

aiii 

BDl 

GII 

BD 

C 

Cl 

Dill 

G1(JV, C)/U(JV) 
G1(AT, R)/0 (N) 
U*(27V)/Sp(A0 
U(p,g)/U(p) x U (q) 
SO(p, q)/SO(p) x SO(g) 
Sp {p, g)/Sp(p) x Sp (q) 
SO(N,C)/SO(N) 
Sp(7V,C)/Sp(iV) 
Sp(A r , R)/U(AT) 
SO*(2A r )/U(A r ) 

U (JV) 

U(N)/0(N) 

U(2AT)/Sp(AT) 

\J{p + q)/V{p) x U (q) 
SO(p + q)/SO(p) x SO (q) 
Sp (p + q)/Sp(p) x Sp (q) 
SO (TV) 

Sp(iV) 

Sp(N)/U(N) 

SO(2N)/U(N) 


Table 1: the large families of symmetric spaces 

The difference from the standard table (37) is that some of the entries of Table 1, namely the spaces of type A, AI and 
All, are not irreducible. They can be made so by dividing out a factor U(l) (R + ) in the compact (resp. noncompact) 
cases. Division by such a factor is analogous to removing the center of mass motion from a mechanical system with 
translational invariance. It turns out that, with a view to superanalytic extensions (cf. Example 2.4), it is preferable 
not to insist on irreducibility but to “retain the center of mass motion”. 

The next subsection introduces super-generalizations of Cartan’s symmetric spaces which have appeared in the 
theory of mesoscopic and disordered single-particle systems and have come to play an important role in that field. 


F. Riemannian symmetric superspaces (definition) 


Let G\ be a complex Lie supergroup that is realized as a group of supermatrices g = ( 900 901 ), with matrix 

elements that take values in a (sufficiently large) parameter Grassmann algebra A = Ao + Ai. If Gc = Gc + Gc is 
the Lie superalgebra of Ga, the Lie algebra of Ga is obtained by taking the even part of the tensor product with A: 
Lie(GA) = A 0 (8> Gc + Ai ® Gc = (A (g) <7c)o- Thus, if {e*, ej} is a homogeneous basis of complex matrices in Gc, an 
element X £ Lie(GA) is expressed by X = z l ei + (Hj with z 1 £ A 0 and (/ £ Ai. 

Let 6 : Ga —» Ga be an involutory automorphism and let Ha C Ga be the subgroup fixed by 6. The decomposition 
into even and odd eigenspaces of 0* : Lie(GA) —> Lie(GA) is written Lie(GA) = Lie(ffA) + Ma- This decomposition 
is orthogonal with respect to the Ad(GA)-invariant quadratic form B : Lie(GA) x Lie(GA) —* Ao, B(AT, Y) := STrXF. 

Both Ga and Ha are supermanifolds with underlying spaces that are Lie groups and are denoted by Gc and He- 
Passing to the coset spaces one obtains a graded commutative algebra H = Hq + Hi of (Grassmann-analytically 
continued) holomorphic sections of the bundle Ga/Ha —► Gc/Hc■ These sections are called (super-)functions 
(on Ga/Ha ) for short. In local complex coordinates z 1 3 * ,..., z p ; ..., ( 9 they are written f(z 1 ,...,z p ;£ 1 ,..., ( q ) = 
Y/ fii—in ( zl > •••! z p )</ 11 •■■ 0 " where the coefficients (z 1 ,..., z p ) take values in A after Grassmann-analytic contin¬ 

uation. For coordinate-independent calculations the alternative notation f(gHA ) or f(g ■ o) is used. In the following 
Gc/Hc is assumed to be connected. 

Every X £ Lie(GA) is associated with a vector field (or even derivation) X : H —> H by 


(■ Xf)(g-o ) 


d_ 

ds 


Ke~ sX g-o). 

0 


( 9 ) 


Here e sX g means the usual product of supermatrices, and the function f{e sX g ■ o) is determined from f(g ■ o ) b; 
Grassmann-analytic continuation. The Lie algebra of even derivations of Ti is a left 7do-module denoted by Dero7d. 


2 We here do not distinguish between the orthogonal groups in even and odd dimension. 

3 Since we always have the option of expanding with respect to the anticommuting parameters that may be containe d in X, 

no information is lost by not considering the full left H-module of supenierivations of 77, cf. the last paragraph of Sec. II B. 
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A notion of supergeometry on Ga/Ha is introduced via a left-invariant tensor field (•, •} : DeroTf x Dero7d —> Hq. 
The details are as follows. Ga acts on Ga/Ha by left translation, : f[g ■ o) i—> f[[hg) ■ o ). The left-translate 
cLTh{X) of a vector field X is defined by the equation Tj*(dTh(X)f) = X(Tj/f), and one requires: 

Th (dT h (X), dT h [Y)) = {X,Y). 

This equation determines (•, •) uniquely within a multiplicative constant. For vector fields of the special form i) one 
obtains 

(X, Y)(g -o) = c o X B ((Adfo)" 1 *)^, (Adfo)"^)^) , 

where the subscript Ma means projection on the odd eigenspace of 0*. Note that since (Ad(g/i) _1 A )= 
Ad(/i) _1 (Ad( 3 ) _1 A')^ v|A for h £ Ha , this is well-defined as a function on Ga/Ha ■ The normalization is fixed 
by choosing cq = 1. 

The metric tensor (•,•) induces a geometry on the ordinary manifold Gc/Hc by restriction (i.e. by setting all 
anticommuting variables equal to zero). Of course, since the groups Gc and He are complex, this geometry is 
never Riemann. However there exist submanifolds in Gc/Hc which are Riemannian symmetric spaces and can be 
constructed by selecting from the tangent space T a [Gc/Hc) a Lie-triple subsystem Ad (i.e. [Ad, [At, Ad]] C At) such 
that the quadratic form B restricted to Ad is of definite sign. It is then not hard to show |l7j that the image of Ad 
under the exponential map A' i—> e x H a is Riemann in the geometry given by restriction of (•, •). Its completion is a 
symmetric space. 

Definition 2.5: A Riemannian symmetric superspace is a pair [Ga/Ha', M) where M is a maximal Riemannian 
submanifold of the base Gc/Hc- 

Remark 2.6: The merit of this definition is that it avoids any use of an adjoint (or “complex conjugation”) of the 
Grassmann variables. ■ 

By the complex structure of Gc/Hc, the tangent space Adc := T 0 (Gc/Dc) decomposes as Adc = Ad + *Ad where 
Ad is taken to be the subspace of Adc on which the quadratic form B is strictly positive. Now observe that, since an 
element g £ Gc is of the form g = diag(g 0 cn 5n)> the group Gc is a Cartesian product of two factors, and the same 
is true for He- Hence, Gc/Hc factors as Gc/Hc = Mg x M/, and Ad is a sum of two spaces: Ad = Ado © Adi, 
which are orthogonal with respect to the quadratic form B. (It may happen, of course, that one of these spaces is 
trivial.) For Z £ Ad, let the corresponding orthogonal decomposition be written Z = X + Y. Then B restricted to 
Ad is evaluated as 

B [Z, Z) = Tr 0 X 2 - TriF 2 , 

where the relative minus sign between traces is due to supersymmetry (STr = Tro — Tiq). The positivity of B on Ad 
is seen to imply X = A't and Y = —Y' (the dagger denotes hermitian conjugation, i.e. transposition in conjunction 
with complex conjugation). 

Given Ga/Ha, the condition that M be Riemann and maximal in Gc/Hc, fixes M uniquely up to two possibilities: 
either T 0 [M ) = Ad, or T a [M) = iA4. In either case, M is a product of two factors, M = Mo x Mi, both of which 
are Riemannian symmetric spaces. In the first case, Mo is of noncompact type and M\ is of compact type, while in 
the second case it is the other way around. We adopt the convention of denoting the compact space by Mp and the 
noncompact one by Mb- 

In view of Cartan’s list of symmetric spaces (Table 1), we arrive at the following table of large families of Riemannian 
symmetric superspaces: 


class 

Ga/Ha 

M b 

Mj? 

A\A 

Gl(m|n) 

A 

A 

Al\AU 

Gl(m 2n)/Osp(m 2n) 

AI 

All 

All\Al 

Gl(m 2n)/Osp(m 2n) 

All 

AI 

aiii aiii 

G1(toi + m 2 ni + n 2 )/Gl(mi|ni) x Gl(m 2 \n 2 ) 

AIII 

AIII 

BD\C 

Osp(w|2n) 

BD 

C 

C\BD 

Osp(w|2n) 

C 

BD 

CII-DIII 

Osp(2m 2n)/Gl(m n) 

Cl 

Dill 

Dill | Cl 

Osp(2m|2n)/Gl(m n) 

Dill 

Cl 

DDI|CII 

Osp(mi + m 2 2m + 2n 2 )/Osp(mi 2m) x Osp(m 2 2n 2 ) 

DDI 

CII 

CII|DDI 

Osp(rm + m 2 2m + 2n 2 )/Osp(mi 2m) x Osp(m 2 2n 2 ) 

CII 

DDI 


Table 2: Riemannian symmetric superspaces 








Although the entries A\A, BD\C and C\BD look extraneous because they are groups rather than coset spaces, they 
fit in the same framework by putting by Ga = G x G and #(< 71 , 52 ) = ( 52 , 51 ), so Ha = diag(G x G) ~ G and 
Ga/Ha -G. 

As far as applications to random-matrix theory and disordered single-particle systems are concerned, the most 
important structure carried by Riemannian symmetric superspaces is their GA-invariant Berezin measure. Such a 
measure always exists by Definition 2.1 and the existence of local coordinates. To describe it in explicit terms, one 
introduces a local coordinate system by the exponential map Ma —> Ga/Ha , Z 1 —> exp(Z)HA■ By straightforward 
generalization (replace the Jacobian by the Berezinian) of a corresponding calculation (cf. |p7| ) for ordinary symmetric 
spaces, one obtains for the principal term of the invariant Berezin measure the expression DZ o J{Z) where DZ = 
dz 1 A ... A d z p ® 8^...8qi denotes the flat Berezin measure on Ma, and if Tz : Ma — ► Ma is the linear operator 
defined by 


f. ad 2n (Z) 

z k ^ +1 ) ! ’ 

the function J(Z) = SDetT^. (Note x 2n /{2n+ 1)! = x~ l sinhx.) A universally valid expression for the anomaly 

in these coordinates is not available at present. 


III. SUPERSYMMETRY APPLIED TO THE GAUSSIAN RANDOM-MATRIX ENSEMBLE OF CLASS C 


The goal of the remainder of this paper will be to demonstrate that Riemannian symmetric superspaces, as defined 
in Sec. IIF, arise in a compelling way when Gaussian ensemble averages of ratios of spectral determinants for random 


matrices are considered in the large-A limit. The example to be discussed in detail will be the Gaussian ensemble 
defined over the symplectic Lie algebra sp(TV), which has recently been identified j24j] as a model for the ergodic limit 
of normal-superconducting mesoscopic systems with broken time-reversal symmetry. 


A. The supersymmetry method: a simple example 


The pedagogical purpose of this first subsection is to illustrate our strategy at a simple example |58| . If u(TV) is the 
Lie algebra of the unitary group in N dimensions, consider on iu(N) (the hermitian N x N matrices) the Gaussian 
probability measure with width v/y/N. Denoting by H the elements of iu(N) and by dH a Euclidean measure, we 
write the Gaussian probability measure in the form dy,(H) = exp(—TVTrU 2 /2v 2 )dH , J d/i(H) = 1. This measure is 
called the Gaussian Unitary Ensemble (GUE) in random-matrix theory. The object of illustration will be the average 
ratio of spectral determinants, 


Z(a,0) 



where a, 0 are complex numbers and a is not in the spectrum of H. Given the generating function Z , the GUE 
average resolvent is obtained by 



z) 1 dy( y H) 


da 


Z(a,0 ) 


cx=(3=z 


We will now show how to compute Z using a formalism that readily generalizes to more complicated situations. 

To avoid the introduction of indices and have a basis-independent formulation, we choose to interpret H as a 
self-adjoint endomorphism H £ End(U) of TV-dimensional complex space V := C N with a hermitian quadratic form 
(x,y) i-> {x,y) v . 

The supersymmetry method starts by introducing “bosonic space” Wb = Wo = C and “fermionic space” Wf = 
Wi = C. Auxiliary space is the Z 2 -graded sum W = Wb ® Wf = C 1 ! 1 . The Cartesian basis of W is denoted 
by e-B = (1,0) and eF = (0,1). Let Hohia (W,V) := Ao (g> Hom(WB,U) + Ai ® Hom(WF,U) where A = Ao + Ai 
is the Grassmann algebra with dimcHom(Wp, V) = N generators. (Grassmann-analytic continuation will not be 
needed here.) Hoim (U, W) is defined similarly, with another Grassmann algebra A. The key idea is to utilize the 
Gaussian Berezin integral over the complex-analytic superspace HomA(W, V) x Horn 3 ^(D, W ). Let D(ip, ip) (with ip £ 
Houla(W, V) and ip £ Horn-^U, W)) denote a translation-invariant holomorphic Berezin measure on this linear space. 
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If 0 b (0b) is the restriction of 0 (0) to a map Wb —> V (resp. V —> Wb), fix a Berezin integral / i—> f E(0, 0) 

by choosing for the domain of integration the subspace M T selected by the linear condition 0 b = 0b (W e adjoint 0g : 
C N —> C being defined by = {z, 0b ■ l)y). Because Houla(W, V) x Hom^(V) W) has complex dimension (2 N, 2TV), 
the integral / E(0,0)/(0,0) does not change its value when / is replaced by the rescaled function / s (0,0) = 
/(s0, s0) (s £ R). Now with End 0 (W) = End(Ws) © End(Wp) and Endi(W) = Hom(WB, Wf) © Hom(WF, Wb), let 

End A (W) := A 0 © End 0 (W) + Ai © Endi(W), 

where A = Ao+Ai is the Grassmann algebra with dimcEndi(W) = 2 generators, and pick A £ End(E),E £ End A (W). 
B corresponds to what is called a 2 x 2 supermatrix in physics. An elementary but useful result is that, if we normalize 
E(0,0) by / E(0,0)exp s 2r Er00^ = 1, the identity 

J D(i/j , 0) exp ^*TryA00 — iSTrwBtfiipj = SDety®vy (A © 1 — 1 © B)~ c (10) 

holds with c = 1 provided that the integral exists. (The parameter c is introduced for later convenience.) When 
A and B are represented by diagonal matrices, verification of © is a simple matter of doing one-dimensional 
Gaussian integrals. The general case follows by the invariance of .0(0,0) under unitary transformations of V and 
“superrotations” in W. 

Now introduce elements Ebb and Eff of Endo(W) by Ebbcb = es, Effcf = eF, and Ebbcf = Effcb = 0. By 
setting A := H and B := ccEbb + /3Eff =: u>, and using 

SDety® )w{H © 1 — 1 © u>) = Det (H — a)/Det (H — f3), 

we get a Gaussian integral representation of Z: 

Z(uj) :=Z(a,/3) = J SDety®ry(E © 1 — 1 © co)~ c d^(H) 

= J E(0,0) J exp ^iTryE00 — iSTrryw00^ d/r(E). (11) 

In the next step, the GUE ensemble average is subjected to the following manipulations: 


exp(zTrE00)d^(E) 

(iTr E00 - NTiH 2 /2v 2 ^ 


exp ' 


dH 


' iu(N ) 


= expTry (00) 2 = exp-^-STrvy(00) 2 
= j DQ exp r*STrQ00 — 7VSTrQ 2 /2u 2 ^ 

JRxiR ^ ' 

=: f D[i(Q) exp(iSTrQ00). 


( 12 ) 


The fourth equality sign decouples the quartic term STrw(00) 2 by introducing an auxiliary integration over Q £ 
End A (W). In order for this Gaussian integral to converge, the integration domain for the BB-part Qbb : Wb —> Wb 
(FF- part Qff : Wf —> Wf) is taken to be the real (resp. imaginary) numbers. By using the relations ( [TT^ . [TT| . [1~2| ) we 
obtain 


Z(u>) = j D(ip,^}) yj exp(jTryE00)d/x(E) J exp — «STiyyw00 

= J J £>(0,0) exp*Try0(Q - w)0 
= J D[i(Q) SDety®jy (ljv © {Q — u>)) c 

= J D^l(Q) SDetw(Q — ui)~ N 

= f DQ exp — NSTi (<Q 2 /2v 2 + ln(Q — w)) . 

JRxiR 


(13) 
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These steps reduce an integral over the N x N matrix H to an integral over the 2x2 supermatrix Q. The large 
parameter N now appears in the exponent of the integrand, so that the Q-integral can be evaluated by a saddle-point 
approximation that becomes exact in the limit N —> oo. By solving the saddle-point equation —Q/v 2 = (Q — w) _1 
and doing an elementary calculation, one obtains Wigner’s semicircle law for the GUE density of states |p9| : 

f Tr :S(E - H)dy(H) = — Wl-(E/2v) 2 , (14) 

J nv 

which will be of use later. 


B. Definition of the Gaussian Ensemble of type C 


Having run through a simple and well-known example, we now treat in detail a less trivial case where the reduction 
to a Q-integral representation requires more care. 

The “physical space” of our model is V = C 2 0 C N . As before, let x i—> x denote complex conjugation, and fix 
a symmetric quadratic form (*,*)y : V x V —> C such that the corresponding hermitian quadratic form ( x,y)v = 
(y,x) v is strictly positive. The transpose and the adjoint of a linear transformation L £ End(H) are defined by 
( x,L T y)v = (Lx,y)v and (x,L'y)v = ( Lx,y)v as usual. 

Consider now the space, P, of self-adjoint Hamiltonians H £ End(H) subject to the linear condition 

H = -CH t C~\ (15) 

where C is skew and C 2 = —1. Clearly, iP is isomorphic to sp(iV) = Cn (the symplectic Lie algebra in 2N dimensions). 
Introducing an orthonormal real basis of V we can represent H by a 2N x 2 N matrix. The explicit form of such a 
matrix is 

b 


H = 


—a 1 


if 


C = 


0 

— 1 jv 


Ijv 

0 


where a ( b) is a complex hermitian (resp. symmetric) N x N matrix. The Gaussian ensemble to be studied is defined 
by the probability measure dy(H) = exp(— NTtH 2 /2v 2 )dH, f dy(H) = 1. For any two A, B £ End(H), 

r 7,2 

/ Tr(AH)Tr(BH) dy(H) = — Tr (AB - ACP T C _1 ) . (16) 

Jixsp(N) 2 AT 

The joint probability density for the eigenvalues of H has been given in (24). 

The physical motivation for considering a Gaussian random-matrix ensemble of the above type (type C) comes 
from the fact |24|] that it describes the ergodic limit of mesoscopic normal-superconducting hybrid systems with 
time-reversal symmetry broken by the presence of a weak magnetic field. To deal with such systems, the Bogoliubov- 
deGennes (BdG) independent-quasiparticle formalism is used. The first factor in the tensor product V = C 2 (g> C N 
accounts for the BdG particle-hole degree of freedom, which is introduced for the purpose of treating the pairing field 
of the superconductor within the formalism of first quantization. The second factor represents the orbital degrees of 
freedom of the electron. H is the Hamiltonian that enters into the BdG-equations, and the relation (|l5[) expresses 
the particle-hole symmetry of the BdG-formalism. 

Our goal is to compute the following ensemble average: 


Z n (ai,..., a„; ft, ...,/?„) = [ J^Det ) 

■ 'isp(N) \H-aiJ 


dy(H). 


(17) 


By the particle-hole symmetry of iJ, Z n is invariant under a reversal of sign for any pair ( cti,Pj ), so no information 
is lost by restricting all on to one half of the complex plane. For definiteness, we require 

Imctj < 0 (* = 1, ...,n). (18) 

All information about the statistical correlations between the eigenvalues of H can be extracted from Z n . For example, 
the probability that, given there is an eigenvalue at E\, there exist n— 1 eigenvalues at E 2 ,..., E n (regardless of the 
positions of all other eigenvalues) is equal to 


Rn(El 


•••, E n )= lim 

E —>0 



n n 

TT 9 

n o 

TT d 

V 7T ) 

}±da 2l 

a 2 i —Ei—ie ^ da 2 l-l 


a 2 i-i=-Ei-ie 


^ ^2n (^ 1 1 • • • 5 ^2n ? -^1 -^1 Z£^...^E'n E n is'). 

The function R n (Ei,E n ) is called the n -level correlation function in random-matrix theory |3S 


(19) 
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C. Symmetries of the auxiliary space 


To transcribe the supersymmetry method of Sec. Ill A to the computation of Z n (which involves n ratios of spectral 
determinants), a simple and natural procedure would be to enlarge the auxiliary space W by taking the tensor product 
with C”. However, on using the formula 


exp(jTr Hipip)dy{H) = exp — - 


f isp(iV) 


(Tr Hipip) 2 dy(H) 


l 




one faces the complication that the second moment f (Tr Hipip) 2 dy(H) then is a sum of two terms, see the right-hand 
side of ©• Consequently, one needs two decoupling supermatrices Q (one for each term). Although this presents 
no difficulty of a principal nature, it does lead to rather complicated notations. An elegant remedy is to modify the 
definition of ip and ip so that ipip shares the symmetry (15) of the BdG-Hamiltonian H. The two terms then combine 
into a single one: 


J (Tr Hipip) 2 dy(H) = ^STr w $iP) 2 , 

which can again be decoupled by a single supermatrix Q. To implement the symmetry ©, we proceed as follows. 

We enlarge the auxiliary space W = Wb ® Wf in some way (left unspecified for the moment) and fix a rule of 
supertransposition Horn \{W,V) —> HoniA(U, W), ip i—> ip T , and Hom^V, W) —> Hom^VF, V), ip i—> ip T . Such a 
rule obeys ip TT = ip a and if rT = crip, where a £ Endo(W) is the operator for superparity, i.e. a(x + y) = x — y for 
x+y £ Wb® Wf = W. It induces a rule of supertransposition EndA(W) —> EndA(W), Q i—» Q T (no separate symbol is 
introduced). Combination with complex conjugation gives a rule of hermitian conjugation f : Endo(W) —> Endo(W). 
Now impose on ip £ Houla(W, V), ip £ Hom^T, W) the linear conditions 


ip = Cip T 7 1 , ip =—-yip T C 1 , 


( 20 ) 


with some invertible even element 7 of Endo(W). The mutual consistency of these equations requires 

7 = 7 T cr. ( 21 ) 

To see that, insert the transpose of the second equation in (|o|) into the first one. Using ip TT = ip a you obtain 
ip = —CC~ lT ipa^ T 'y~ 1 . Since CC~ lT = —1 and <T 7 T = 7 T ct, Eq. ( pl| ) follows. The consistency condition can be 
implemented by taking Wb = Wf = C 2 ® C", see below. By multiplying the equations (^0|) we obtain 

1 pip — — C(t/>'i/>) T C -1 , ipip = —'y{ipip) Tr t~ 1 - ( 22 ) 


The first equation is the desired symmetry relation allowing us to combine terms. To appreciate the consequences of 
the second equation, note that by the fourth step in (El the symmetries of ipip get transferred onto Q, so that the 
latter is subject to 

Q = -7Q T 7 _1 - (23) 

This symmetry reflects that of the BdG-Hamiltonian H, see (ED- The linear space EndA(W), when given a Lie 
bracket by the commutator, can be identified with gl(2n|2n) = Lie(Gl(2n|2n)). As 7 is supersymmetric (7 = 7 T cr), 
( p3| ) fixes an osp( 2 n| 2 n)-subalgebra. 

7 is not unique. For definiteness we choose it as follows. Let {Eij }aj=i,...,m be a canonical basis of End(C M ) 
satisfying EijEki = SjkEu (here M = 2 or M = n). For M = 2 define the Pauli spin operators a x = E 12 + E 21 , 
o y = —iEr 2 + 1 E 21 , and a z = f?n — The usual rule of supertransposition on EndA(W) is given by (y, v = 1,2 
and i,j = 1 ,..., n) 


( Ebb <8 E^ v ® Eij) T = Ebb ® E vll ® Eji, 
(Ebb ® E^ v ® Eij) T = Ebf ® E VI1 0 Eji, 
With these conventions, one possible choice for 7 is 

7 = -Ebb 0 7b + -Eff <8 7f where 


(Ebf 8 E M1/ 0 Ey) T = — Efb 8 E v/X 0 Eji, 
(Eff 8 E M „ 0 Ey) T = Eff 8 E vtl 0 Eji. 

7B = 0 In, 7F = i<y y 8 In- 


(24) 


This is the choice we make. 
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D. Gaussian Berezin integral 


To repeat the steps of Sec. Ill A and derive a Q-integral respresentation for the generating function Z n , we must 
first generalize the basic identity (|l(i|), whose left-hand side is 


D(ip, if) exp [iAny Aifif — iSTi'n/i? 


By (22) we have 


TrAifip = Tr (if if ) T A T = ^Tr (A — CA T C ^ifif, 
STr Bifift = STr('i fif) T B T = |Tr (B — r yB T 'y~ 1 )ifif. 

In view of this we demand that A and B satisfy: 

A = ~CA T C-\ B = 


(25) 


(26) 


When carrying out the calculation HD we need to apply the identity (jn]) twice, the first time with A = H, B = u>, 
and the second time with A = 0, B = u> — Q. In order for (^ 6 |) to be satisfied with these identifications, we choose to 
set 


n n 

CO = Ebb ® cr z ® ^ ctiEu + E FF <g> a z <g> ^ PjEjj. 

i= 1 3 =1 


The presence of the factor a z = diag(+l, —1) reverses the sign of the a* and ,6j on that subspace where a z acts by 
multiplication with minus one. As the imaginary parts of the at control the convergence of the integral, this sign 
reversal affects the correct choice of integration domain for i/jb and ipB- To ensure convergence of the integral (p5|), we 
require ImSTrwf/'f/; < 0. This inequality is achieved by imposing the condition -0 b = (o z g 1 n)^, which is compatible 
with C = i(Jy g Iat, i/jb = Ci)bTb > an d 7b = &x g In- 

Lemma 3.1: Let D(if,if) denote a translation-invariant holomorphic Berezin measure on the subspace of 
Hom^IT, V) x Houia(F, W) defined by @). Fix the integration domain by 0 b = (c z g l n )V4> an d normalize 
Zi>(0,0) so that / D(if, 0) exp(— s 2 Trifif) = 1 (s el). Then if A £ End(F) and B £ EndA(IF) are diagonalizable 
and satisfy the linear conditions (^6|), the identity ([To|) holds with c = 1/2 provided that the integral exists. 

Proof. Assume that A and B are represented by diagonal matrices 


N 


A = a z gJ2 x ^ B = E B B®* z ®Y. z i E n+ E ™®°*®Y. Vj E jj > 

3 =1 3=1 


i=1 


which conforms with (26). The right-hand side of (0) then reduces to 


N n 


SDet V „(A ® 1 - 1 * = n II . 

1= 1 7=1 ( Xi ~ z i)( Xi + z i> 


(27) 


To evaluate the left-hand side write 


IpB = 


V'F = 


a P 

7 6J> 


where a, 6 , c, d (a,/3, 7 , <5) are complex N x n matrices with commuting (resp. anticommuting) matrix elements. The 
constraint if = — 7 ^ T C _1 results in 


if B = 


-d T b T 
-c T a T 


ifF = 


-S T f3 T 
7 t — a T 


and the reality condition ?/>b = (07 g l n )')/| 3 means d = —a and c = b. The exponent of the integrand is expressed by 
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N n 

— |Tr B^i/j = E £((**- Z j )V j ij V'ij ( J- t + 7/ ) b'i j ^i j 

i=1 j =1 

T (*U yj)Pijlij) ■ 

Doing the Gaussian integrals one gets a result that is identical to (|27]), which proves the Lemma for diagonal A and 
B. The general case follows by the invariance properties of D{ip,ip). 

Remark: The condition of diagonalizability can of course be weakened but we won’t need that here. ■ 

To apply Lemma 3.1 to our problem, note 


SDety^jy {H 0 1 


n 

1 0 w) 1/2 = Dety 
2=1 


f {H - ai)(H + a t ) 


1/2 


\jDet 

2=1 


(H-aA 

\H-h) 


where in the second step we used the invariance of the ratio of determinants under H i—> —H, which is due to the 
particle-hole symmetry H = —CfL T C _1 . Moreover, note 

SDety®vv (1 0 (Q — lo)) 1 ^ 2 = SDetw(Q — w) N ■ 

The previous calculation ( 00 ) thus formally goes through with c = 1/2, and isp(TV) for iu(N), and we arrive at 
the following representation of the generating function: 


Z n {u)) 


DQ exp —NSTr (Q 2 /2v 2 + ln(Q — w)) , 


(28) 


where the supermatrix Q = ^ ^ BB ^ BF j is subject to (0). To make this rigorous, we have to specify the integration 
domain for Q and show that the interchange of the {ip, ip)- and Q-integrations is permitted. 


E. Choice of integration domain 

If the steps ( |TT| - [Uj| ) are to be valid, we must arrange for all integrals to be convergent, at least. This is easily 
achieved for Qff, the FF-component of Q, but requires substantial labor for ipB, ipB and Qbb- Consider Qff first. 
Since —STrQ 2 = — TrQg B + Tr+ nilpotents, we want TtQffQff < 0, which leads us to require that Qff be 
antihermitian. Combination with (|23| ) gives 

Q ff = 7fQff7f = Qff 

where 7 f = icr y 0 1„, see ([24]). The solution space of these equations is sp(n), the symplectic Lie algebra in 2 n 
dimensions. Thus we choose U := sp(n) for the integration domain of QfFj and of course the integration measure is 
taken to be the flat one. 

The choice of integration domain for Qbb is a much more delicate matter and will occupy us for the remainder of 
this section. Recall, first of all, that the convergence of 

J D{ip, ip) exp ^ iTrHipip — iSTrunpip^j 


requires taking ipB = AP'n where /3 := cr z 0 l n cancels the minus signs that multiply the imaginary parts of the 
parameters a, in u. To ensure the convergence of 

J D{ip, ip) exp iTri/>(Q — Lo)ip, 

one is tempted to choose Qbb in such a way that ReTY ipQip = 0. Unfortunately, when this condition is adopted one 
gets Qbb = PQbbPi which causes TrQ| B = ^QbbPQbbP t° be of indefinite sign, so that the integral over Q does 
not exist. 

A way out of this difficulty was first described by Schafer and Wegner ( 2 ) in a related context. We are now going 
to formulate their prescription in a language that anticipates the geometric structure emerging in the large -N limit. 
To simplify the notation, we put Qbb = iZ. What we need to do is investigate the expression 
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(29) 


exp(—A r TrQg B /2t! 2 + «Tt(3bbV’bV’b) = exp(7VTrZ 2 /2i; 2 - TtZ^b^b)- 
The conditions on Qbb translate into 

Z = -Tb^Tb 1 = -P&p- 1 . 

Because 7 b = <r x ® l n is symmetric, the solution space of the first equation is a complex Lie algebra Gc — so(2n, C). 
The matrix representation of an element Z G G c is of the form ^ ^ ^ where B and C are skew. The second 

equation (Z = —/3 iU/3 _1 ) means A = —A^ and C = B\ which fixes a real form Q = so*(2 n) of Gc = so(2n,C). This 
real form is noncompact (i.e. G = Lie(G) with G a noncompact Lie group), which is what causes all the trouble and 
is forcing us to work hard. Its maximal compact subalgebra 1C is the set of solutions of X = /3X(3~ 1 in G■ From 

X = ^ ^ and A = —A^ we see that K. ~ u(n). 

To display clearly the general nature of the following construction, we introduce a symmetric quadratic form 
B : Gc x Gc ► C by B(X, Y) = TrA'F. The Cartan (orthogonal) decomposition of G with respect to this quadratic 
form is written G = 1C ® M. An element Y of M satisfies Y = —/3Yf 3 _1 . From this in conjunction with the equation 
fixing 1C (X = +PX/3~ 1 ) one deduces the commutation relations 

[M,M\ C/C, [K,M\CM, [/C,/C]c/C. (30) 

Note that the elements of M are hermitian while those of K, are antihermitian. We will also encounter the complexified 
spaces K-c = 1C + ilC and Me = M + iM. They, too, are orthogonal with respect to B and satisfy the commutation 
relations @). The element (3 = a z l n satisfies f3 = —Tb^Tb 1 and can therefore be regarded as an element of Gc- 
Moreover, (3 £ ilC C Gc- 

Now we embed G = 1C © M into Gc by a map fo, 

fib ■ 1C x AI—> Gc, 

(X , Y)» fi b {X , Y) = b x (X + e Y Pe~ Y ), 

where b ^ 0 is some constant that will be specified later. 

Lemma 3.2: fib(lC x M) is an analytic manifold without boundary, and is diffeomorphic to G■ 

Proof. Analyticity is clear. To prove the other properties, we first establish that fi b is injective. For that purpose, 
we write ef /3e~' 1 = e ad ^ F ^/3 where ad {Y)f3 = [Y,/3] is the adjoint action on Gc- Decomposing the exponential function 
according to exp = cosh + sinh, we write fi b = fi+ + fi~ where 

fi+(X,Y) = bx (A' + coshad(Y)/3), 
fi-(X,Y) = bx sinlrad(Y)/3. 

From the commutation relations (^) and (3 G iK. we see that fi± takes values fi+(X, Y) G K-c and fi~{X, Y) G Me- 
Since Gc = © Me (direct sum), injectivity is equivalent to the regularity of the maps X i—> fi + (X,Y) (with Y 

viewed as a parameter) and Y i—> fi_[X,Y). The function fi + ( A, •) = X + const is obviously regular. By Y = Y' 
the element Y is diagonalizable with real eigenvalues. The regularity of fi- then follows from sinh : R —» R being 
monotonic and Y i—> ad (Y)/3 being regular. This completes the proof that fi b is injective. The injectivity of fi b means 
that fi b (!C x M) is diffeomorphic to G = 1C © M. This in turn means that, since G has no boundary, fi b (!C x M) has 
no boundary either. ■ 

We are now going to demonstrate that fi b (lC x M) for any b > 0 may serve as a mathematically satisfactory 
domain of integration for the variable Z in (|29|). We begin by investigating the quadratic form TrZ 2 = B (Z,Z) 
on fi b (lC x M). For this we set Z = Z + + Z_ with Z± = fi±(X,Y). Using B (Z + ,Z~) = 0 (recall ICc-LMc)-, 
B(ad(F)A, B) = — B(A, a,d{Y)B) and cosh 2 — sinh 2 = 1, we obtain 

B(Z,Z)/b 2 = B(A, X) + 2B(A,coshad(F)/3) +B(/3,/3). 

The antihermiticity of X G 1C gives B(A', X) < 0. In contrast, coshad(F)/3 G ilC is hermitian, so B(A, coshad(U)/3) G 
It follows that exp(A r TrZ 2 /2u 2 ) = exp(NTrfi b (X,Y) 2 /2v 2 ) is decaying with respect to X and oscillatory w.r.t. 
Y. 

We have not yet made any use of b > 0 yet. This inequality comes into play when the coupling term 
-TrZfiBfiB = ~B(Z, fisfis) = —bB(X, fisfiB) - bB(e y f3e~ y ,^b'0b) 
is considered. From (^2|) and = (3-fi^ we see that fiBtpB satisfies 
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= -7bW’b'0b) T 7b 1 = +/3(iPb4’b)^P \ 

so V’b^b € it/. Since B is real-valued on Q x Q, the term B(A', ^b^b) is purely imaginary. The other term, 

-M3(e 5 /3e _1 ,'tpBipB,) = -6Tr((/> B e 2y 4) < 0 

is never positive if b > 0. Hence the real part of the exponential in ( |29| ) is negative semidefinite for Q = iZ £ i(fb(K,x M.) 
and b > 0. As a result, the integrals over Q and if,if converge if the integration domain for Q is taken to be 
i(f>b(IC xM)xU (b > 0). Because i<j>b(IC x Ad) x U is an analytic manifold without boundary and Cauchy’s theorem 
applies, we may perform the shift of integration variables that is implied by the fourth equality sign in & . Moreover, 
the presence of the nonvanishing imaginary parts of of the parameters cq in u> ensures uniform convergence of the 
(if, if)- integral with respect to Q, so that we may interchange the order of integration (the second equality sign in 
(13)). And finally, any breakdown of diagonalizability of Q — uj occurs on a set of measure zero, so that the identity 
( lCf) (Lemma 3.1) may be used, and all steps leading to ( |28| ) are rigorous. In summary, we have proved the following 
result. 

Theorem 3.3: For V = C 2 (g) C N and W = C 1 ! 1 ® C 2 (g> C" define the generating function 

Z n ,N(w) = f SDety 0 vv (H <8 1 — 1 <g> w) -1 ^ 2 exp (— NTrH 2 /2v 2 ) dH, 

J ZXsp(iV) 

n n 

tv = Ebb <8> cr z <8 ^ aiE u + E FF ® cr z <8 ^ PjEjj (Ima, < 0). 

*=i i=i 

Let DQ denote a translation-invariant holomorphic Berezin measure of the complex-analytic superspace osp(2n|2n). 
Then for all N £ N, n £ N and b > 0, DQ can be normalized so that 

Z n ,N(w) = [ DQexp—NSTr (Q 2 /2v 2 + ln(Q — a/)) , (31) 

J i<p b (/CxM)xU 

where U = sp(n), K, ~ u(n), M. is determined by /C © M = so*(2n), and <fb(X, Y) = b (X + Ad(e r )(tr z ® l n )). ■ 

We conclude this subsection with a comment. In the literature a parameterization of the form Q = TPT -1 (cf. [^0[) 
has been very popular. In our language, this factorization amounts to choosing for the integration domain of Qbb 
the image oi ip : Q = 1C (B Xi —> G, X + Y e Y Xe~ Y . This is not a valid choice as <p(G) does have a boundary, 
namely the light cone {Z\B(Z, Z) = 0} in Q , so that shifting of integration variables is not permitted. (However, it 
turns out that the error made becomes negligible in the limit N oo, so that the final results remain valid if that 
limit is assumed.) 


F. Saddle-point super manifold 

The result (|3^) holds for all N £ N. We are now going to use the method of steepest descent to show that in the 
limit N —> oo, the integral on the right-hand side reduces to an integral over a Riemannian symmetric superspace of 
type Hill | CL 

With our choice of normalization, the mean spacing between the eigenvalues of H scales as Af -1 for N —> oo, 
see ©■ We are most interested in the eigenvalues close to zero as their statistical properties describe those of the 
low-lying Bogoliubov independent-quasiparticle energy levels of mesoscopic normal-superconducting systems [Q. To 
probe their statistical behavior, what we need to do is keep w = Nu/nv (i.e. u> scaled by the mean level spacing) 
fixed as N goes to infinity. In this limit to ~ 0(1/N) can be treated as a small perturbation and we may expand 
AfSTrln((3 — uS) = ASTrIn Q — TruSTrQ" 1 :!) + 0(1/N) if Q exists. 

To evaluate the integral ([ll]) by the method of steepest descent, we first look for the critical points of the function 
NF(Q) = NSTr(Q 2 /2v 2 + InQ). These are the solutions of 

F'(Q) = Q/v 2 + Q~ l = 0, 

or Q 2 = — v 2 . The solution spaces, the so-called “saddle-point supermanifolds”, are nonlinear subspaces of osp(2n|2n), 
which can be distinguished by the eigenvalues of Q. Of these supermanifolds, which are the ones to select for the 
steepest-descent evaluation of the integral (D? 

To tackle this question, we start out by setting all Grassmann variables to zero. The BB-part of the saddle-point 
manifold(s) is uniquely determined by the forced choice of integration domain itfb(IC x M.) and by analyticity. This is 


16 



because the saddle-point manifold must be deformable (using Cauchy’s theorem) into the integration domain without 
crossing any of the singularities of SDet(Q — u)~ N ', and by inspection one finds that this condition rules out all 
saddle-point manifolds except for one, which is i<f > v (0 x Ad), the subspace of the integration domain icf>b()C x A4)\b= v 
obtained by dropping from Q = K, ® M. the K, degrees of freedom (these are the directions of steepest descent). By 
an argument given in the proof of Lemma 3.2 we know that i<p v (0 x M) is diffeomorphic to M. On general grounds 
the latter is diffeomorphic to a coset space G/K by the exponential map A4 —> G/K, Y i—> e Y K ; where in the 
present case G = {g £ Gl(2n,C)|.g = 7Bff” lT 7B 1 = /3g -lt /3 -1 }, and K = {k £ G\k = /3kf3~ 1 } (on setting g = exp Z, 
k = exp A and linearizing, we recover the conditions Z = — 7 B^ T 7 g 1 = ~PZ^P~ 3 defining Q and the condition 
X = pXp- 1 fixing the subalgebra 1C). We already know Q = so*(2 n) and K. — u(n), so G = expt/ = SO*(2n) and 
K = exp /C = U(n). Because K is a maximal compact subgroup, the coset space G/K is a Riemannian symmetric 
space of noncompact type. In Cartan’s notation, G/K = SO*(2n)/U(n) is called type DHL For better distinction 
from its FF-analog, we will henceforth denote G/K by G/K-q. 

We turn to the FF-sector. Since SDet(Q — to)~ N does not have poles but only has zeroes as a function of Q ff, 
analyticity provides no criterion for selecting any specific solution space of the saddle-point equation Q p F = — v 2 . 
Instead, the determining agent now is the limit N —> oo. From ( |3l|) it is seen that integration over the Gaussian 
fluctuations around the saddle-point manifold produces one factor of iW 1 (N +1 ) for every commuting (resp. anti¬ 
commuting) direction of steepest descent. Therefore, the limit N —> oo is dominated by that saddle-point manifold 
which has the minimal transverse (super-)dimension dg — dp. A little thought shows that the transverse dimension 
is minimized by choosing Qff to possess n eigenvalues +iv and n eigenvalues —iv. Thus, the dominant saddle-point 
manifold is unique and contains the special point qo := ivP (P = a z (g> 1„ now acts in the fermionic subspace). 

Recall that the integration domain for Qff is a compact Lie algebra U = sp(n). The corresponding Lie group 
U = Sp(ra) operates on U by the adjoint action Ad(u) : U —> U, X uXu -1 . Because the saddle-point equation 
Qff = ~^ 2 Qff invariant under this action, the FF-part of the (dominant) saddle-point manifold can be viewed as 
the orbit of the action of Ad(G) on the special point qo £ U. Let Rp be the stability group of qo, i.e. Kp = {k £ 
Ujkqok -1 = qo}- By Ad(Kp)q 0 = qo the orbit Ad(G)go is diffeomorphic to the coset space U/Kp. Arguing in the 
same way as for the BB-sector, one shows that Kp ~ A'b — U(n). Hence U/Kp = Sp(n)/U(n), which in Cartan’s 
notation is a compact Riemannian symmetric space of type Cl. 

We are finally in a position to construct the full saddle-point sitpermanifold. Recall, first of all, that Q is subject to 
the condition Q = — r yQ Tr y~ 1 , which defines an orthosymplectic complex Lie algebra Q\ '■= osp(2n|2?r) in EndA(W). 
The solution spaces in Q\ of the equation Q/v 2 + Q^ 1 = 0 are complex-analytic supermanifolds that are invariant 
under the adjoint action of the complex Lie supergroup G\ := Osp(2n|2n). They can be regarded as Ad(GA)-orbits 
of elements Qo £ Lie(GA) that are solutions of (Qo) 2 = —v 2 . From the above analysis of the BB- and FF-sectors, we 
know that the saddle-point supernranifold that dominates in the large -N limit is obtained by setting Qo = ivY z where 
= 1 b|f ® P = (Abb + Aff) ® <g> 1„. If Ha is the stability group of Qo, the orbit Ad(GA)Qo is diffeomorphic 

to the coset space Ga/Ha- From yE z + = 0 and the equation hE 2 /i -1 = E z (or, equivalently, h = E z hE 2 ) for 

h £ Ha one infers Ha — Gl(n|n). Hence the unique complex-analytic saddle-point supermanifold that dominates the 
large-AT limit is Ga/Ha — Osp(2n|2n)/Gl(n|n). 

Turning to the integral (31) we note the relations STrQp = —w 2 STrl = 0 and InSDetQo = lnl = 0. These imply 
that the function F(Q) = STr(Q 2 /2v 2 + ln Q) vanishes identically on Ad(GA)Qo- Hence the exponent of the integrand 
in (31) restricted to Ga/Ha is 


7rr>STrQ 1 w 


Ga/H a 


+ 0(1/N) = -ArB (w, Ad( 5 )E 2 ) + 0{1/N). 


To complete the steepest-descent evaluation of © we need to Taylor-expand the exponent of the integrand up to 
second order and do a Gaussian integral. By the Ad (Ga) -invariance of the function NF(Q) it is sufficient to do this 
calculation for one element of the saddle-point supernranifold, say Q = Qq. Putting Q = Qo + Z (Z £ Qa) we get 

N 


NF{Qo +Z) = ^STr (Z 2 + ZY z ZY z ) + 0(Z 3 ). 
2v z 


Now we make the orthogonal decomposition Qa = Lie(AA) + A^a, Z = X + Y, where Y = —E 2 FE 2 are the degrees 
of freedom tangent to the saddle-point supernranifold, and X = +S 2 XE 2 are the degrees of freedom transverse to 
it. The translation-invariant Berezin measure DZ of Qa factors as DZ = DYDX. We thus obtain the transverse 
Gaussian integral 


J DX exp (~NSTrX 2 /v 2 + O(N 0 )) . 


The integration domain for X is HCb x /Cf — 'iu(n) x u(?z). By dim Lie(AA) = (P, q) and p = q, this integral reduces 
to a constant independent of N in the limit N —> oo. 
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What remains is an integral over the saddle-point supermanifold itself. Since DY is the local expression of the 
invariant Berezin measure of Ga/Ha at Ad(e^ )Qo|y=o = Q o we arrive at the following result. 

Theorem 3.4: If Dgn is a suitably normalized invariant holomorphic Berezin measure of the complex-analytic 
supermanifold Ga/Ha ~ Osp(2n|2n)/Gl(n|n), 

lim Z nt Ar{^vu>/N) = / Dgnexp — Ad(g)E 2 ) (32) 

N ^°° ’ JMb X Mf 

where E z = 1 B |f ® 07 ® l n , M B — SO*(2n)/U(n), and Mp ~ Sp(n)/U(n). 

Remark 3.5: This result expresses the generating function for N —> oo as an integral over a Riemannian symmetric 
super space of type DIII|CI (see Tables 1 and 2) with m = n. 

In Q] the n-level correlation function R n is calculated exactly from (|3^) for all n. 


IV. OTHER SYMMETRY CLASSES 

There exist 10 known universality classes of ergodic disordered single-particle systems. These are the three classic 
Wigner-Dyson classes (GOE, GUE, GSE), the three “chiral” ones describing a Dirac particle in a random gauge field 
(chGUE, chGOE, chGSE), and the four classes that can be realized in mesoscopic normal-superconducting (NS) hybrid 
systems. In Ref. [||j] it was noted that there exists a one-to-one correspondence between these universality classes and 
the large families of symmetric spaces (with the exception of the orthogonal group in odd dimensions). Specifically, 
the Gaussian random-matrix ensemble over the tangent space of the symmetric space describes the corresponding 
universality class, in the limit N —> oo. In the notation of Table 1 the correspondences are A <-> GUE, AI <-> GOE, 
All <-> GSE, AIII <-> chGUE, BDl <-> chGOE, Oil <-> chGSE, and the four NS-classes correspond to C, D , Cl, and 
Dill. 

We have shown in detail how to use the supersymmetry method for the Gaussian ensemble over Cjv = sp(IV), the 
tangent space of the symplectic Lie group. There are nine more ensembles to study. We will now briefly run through 
all these cases, giving only a summary of the essential changes. 


A. Class D 


Recall the definitions given at the beginning of Sec. Ill B| and replace the symplectic unit by C = cr x ® ljv- What 
you get is a Gaussian random-matrix ensemble over Dn = so(2 N), the orthogonal Lie algebra in 2 N dimensions. The 
explicit form of the Hamiltonian is 


H = 




where a ( 6 ) is complex hermitian (resp. skew). The treatment of this ensemble closely parallels that of type C. 
A change first occurs in the consistency condition for 7 , which now reads 7 = — 7 T cr (instead of 7 = + 7 T cr) by 
CC~ l = +1. The extra minus sign can be accommodated by simply exchanging the BB- and FF-sectors (73 <-> 75 -). 
The linear constraint Q = ^ 7 Q T 7 -1 again defines an osp(2n|2?r) Lie algebra, the only difference being that the 
BB-sector is now “symplectic” while the FF-sector has turned “orthogonal”. Everything else goes through as before 
and we arrive at the statement of Theorem 3.3 with U ~ so(2n), K. — u(n), and K, © M ~ sp(n, R). 

A novel feature arises in the large-IV limit, where instead of one dominant saddle-point supermanifold there now 
emerge two. One of them is the orbit with respect to the adjoint action of Osp(2?r|2n) on Qq = M'1 B | F ® cr z ® l n as 
before, and the other one is the orbit of 


( n 

Q 1 = ivE B B ®cr z <g> 1„ + ivE-pF <8> o~ <E> I Eu - ^ Eu 

V »= 2 

(The orbits of Qo and Q\ are disconnected because the Weyl group of so(2n) is “too small”.) Consequently, the right- 
hand side of Theorem 3.4 is replaced by a sum of two terms, one for each of the two saddle-point supermanifolds. 
The integral is over a Riemannian symmetric superspace of type CIl-DIII (m = n) in both cases. 
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B. Class Cl 


Let V = C 2 ®C JV carry a hermitian inner product (as always), and consider the space, P, of self-adjoint Hamiltonians 
H £ End(V) of the form 

H = H t = -CHT'C- 1 = (j) where C = ia y <g> 1 N = ^ ^ . 

The N x N matrices a and b are real symmetric. It is easy to see |25|] that P is isomorphic to the tangent space of 
the symmetric space Sp(iV)/U(iV) (type Cl). A Gaussian measure dg{H) on P is completely specified by its first two 
moments, f p Tr (AH)d/j,(H) = 0 and 

r ,,2 

/ Tr(AH)Tr(BH)dg(H) = — Tr (A(B + B T ) - AC(B + B T )C~ 1 ) . 

J p 4 N 

To deal with the random-matrix ensemble defined by this measure, we take W = C 1 ! 1 <S> <C 2 <S> C 2 ® C n . Recall 
ip £ Hom>(W, V) and ip £ Hom^ (V, W). The symmetries of H are copied to ipip by imposing the linear conditions 

ip = Cip T 7 _1 , ip = —'yip T C~ 1 -, ip = tjj T T~ 1 , ip = mp T . 

In order for these conditions to be mutually consistent, r, 7 £ End 0 (IT) must satisfy 

t t -1 -1 

7 = 7 < 7 , r = r cr, 7 r = — T 7 

Without loss, we take 7 and r to be orthogonal. The consistency conditions can then be written in the form 

7 = a = t 2 , 7 T + r 7 = 0 . 

If G1(W) ~ Gl(4n|4n) is the Lie supergroup of regular elements in Endy(W), the equation 7 s = a in conjunction 
with g TT = aga means that the automorphism 7 : G1(W) —> G^W) defined by 7 [g) = 75 _lT 7 _1 is involutory. The 
same is true for f defined by f(g) = rg _lT r _1 and, moreover, 7 and f commute by yr + ry = 0. For definiteness we 
take 


7 = E BB ® 7b + Eff ® 7f, 7b = cr x ®a z ® 1„, 7f = icr y (g> 1 2 ® 1„, 

t = E bb ® t b + E ff ® Tp, t b = 1 2 ® a x <g> 1„, r F = cr- (g) ia y ® l n . 

(This choice is consistent with ip B = /3ip B , (3 = a z <g> I 2 <g> l n .) Let 

Q := {Q £ Endy(IT)|<2 = -yQ T 7 _1 = +t<2 T t -1 } 

be the subspace distinguished by the symmetry properties of ipip. The group G1(W) acts on Q by Q 1 —> gQg~ x . We 
now ask what is the subgroup Gy of G1(W) that leaves the symmetries of Q invariant (the normalizer of Q in Gl(IT)). 
Lemma 4-1: Gy is isomorphic to Osp(2n|2n) x Osp(2n|2n). 

Proof. The conditions on g £ Gy can be phrased as follows: 

7 = 91 g T , T = 9 t 9 T - 

Equivalently, Gy can be described as the simultaneous “fixed point sef’Qof the involutory automorphisms 7 and t. 
We first describe the fixed point set of 7 o t , which acts by (7 o t)( g) = ege~ l where e = —iyr -1 . From the explicit 
expression e = 1b|f < 8 > cr x (g> a y ® l n we see that e has 4n eigenvalues equal to +1, 4n eigenvalues equal to —1, and 
these are equally distributed over the bosonic and fermionic subspaces. Hence the subgroup of G1(W) fixed by 70 -f is 
isomorphic to G + x G_ where G + ~ Gl(2n|2n) ~ G_. Denote the embedding G+ x G_ —> G1(W) by ip(g+, g~) = g. 


4 For lack of a better word we borrow the terminology from manifolds. Of course what is meant here is the supermanifold of 
solutions in G1(W) of the nonlinear equations g = 7 (g) = f (g). 
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The group Ga is the fixed point set of f (or, equivalently, of 7 ) in y?(G+ x G_) (f commutes with 70 f and therefore 
takes <p(G+ x G_) into itself). Note er = —re, e _lT = —e, and for g G y>(G+ x G_) do the following little calculation: 

sr(g) = erg~ lT t~ x = -reg~ lT r _1 = r(eg)" lT r ^ 1 = f (eg). 


Combining this with ey;(g + ,g_) = <^(g + ,—g_) one infers that t acting on y>(G+ x G_) is of the form f</?(g+,g_) = 
<p(r + (g + ), r_(g_)). By a short calculation (work in an eigenbasis of e) one sees that the involutory automorphisms 
Tj : Gl(2n|2n) —> Gl(2n|2n) (i = ±) are expressed by fj(g) = ng -1 r ” 1 with supersymmetric n (t* = T l T cr). It follows 
that fj fixes an orthosymplectic subgroup of Gi ~ Gl(2n|2n), so Ga ~ Osp(2n|2n) x Osp(2n|2n) as claimed. 

Corollary f.2: The space Q is isomorphic to the complement of osp(2n|2?z) ® osp(2n|2n) in osp(4n|4?z). 

Proof. The solution space in EndA(IF) of Q = — 7 Q T 7 _1 is an osp(4n|4n) algebra. Implementing the second condi¬ 
tion Q = +tQ t t ~ 1 amounts to removing from osp(4n|4?z) the subalgebra fixed by X = By linearization 

of the conditions g = 7 (g) = f(g), this subalgebra is identified as Lie(GA) — osp( 2 n| 2 n) ® osp( 2 n| 2 n). ■ 

The Gaussian integral identity © continues to hold, albeit with a different value of c = 1/4. The proof is 
essentially the same as before. 

Since Q is not a Lie algebra, the description of the correct choice of integration domain for the auxiliary variable 
Q is more complicated than before. In the FF-sector we take U := {Qff G QffIQff = —Qf F }- By Corollary 4.2, 
sp(2n) ~ (sp(n) ® sp(n)) ®U. To deal with the BB-sector we introduce the spaces 

g = {X£ gl(2n,C)|X = = -tbA '^ 1 = -/^/T 1 }, 

M = {Y G G\Y = -pYp- 1 }, V ± = {X G Qbb|A = -/JAt/T 1 = ±(3X(3~ 1 }. 


where /? = cr z ® I 2 < 8 > In- The Lie algebra Q is a noncompact real form of the BB-part of Lie(GA). By (3 G iV + and 
the commutation relations [M.,V + ] C P~ and [M,V~] C V + , we have an embedding 


<t>b '■ P + x A4—> Qbb = P^ + Pc > 
{X,Y)^bx (x + e ad(Y) /3) . 


Similar considerations as in Sec. 1IIE show that all integrals are rendered convergent by the choice of integration 
domain 4>bi/P + x A4) xU (b > 0) for Q. With this choice we again arrive at Theorem 3.3. 

The large-TV limit is dominated by a single saddle-point supermanifold, which can be taken as the orbit of Q 0 = ivY z 
(here E z = 1 


B|F 


I 2 < 8 > l n ) under the adjoint action of Ga- This orbit is diffeomorphic to Ga/Ha where 
= E z }. The stability group Ha can equivalently be described as the fixed point set of 


a z ® 

H\ = {h G Ga|/iE z /i _1 

E z : Ga —► Ga, E-(g) = E~gE z . By the relations E z = El = —yE z 7 _1 = tE z t _1 (E z G Q), the element E z 
anticommutes with e = —*yr _1 , and E z commutes with 7 o f. These relations are compatible with the existence of 
an embedding <j) : Osp(2n|2n) x Osp(2n|2?r) —> Gl(kF) such that (E z o </>)(g + ,g_) = </>(g_,g + ). (Such an embedding 
is easily constructed.) Hence Ha — dia,g(Osp(2n|2n) x Osp(2n|2n)) ~ Osp(2n|2n). In this way we arrive at Theorem 
3.4 with Ga/Ha — Osp( 2 n| 2 n), and the maximal Riemannian submanifold Mb x Mp where Mb — SO(2n, C)/SO(2n) 
and Mp ~ Sp(n) (type P|G).^j 


C. Class Dill 


Consider for V = C 2 (g> C 2 ® C N the linear space 

P = {H G End(V)| H = H^ = -C£T t C _1 = +TH t T~ 1 }, 

where C = a x ® I 2 ® In and T = I 2 ® icr y ® 1 n- It has been shown j2^] that P is isomorphic to the tangent space of 
SO(47V)/U(27V) (a symmetric space of type Dill). Introducing an orthonormal real basis of V, we can represent H 
by a 47V x 47V matrix. If C and T are given by 


5 Since the orthogonal group here always appears with an even dimension, we use the simplified notation D\C, instead of 
BD\C as in Table 2. 
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/ 0 

0 

1 N 

0 ^ 


( 

0 

ljv 

0 

° \ 

0 

0 

0 

ljv 

, T = 


-ljv 

0 

0 

° 

1 N 

0 

0 

0 


0 

0 

0 

ljv 1 

V 0 

1 N 

0 

0 ) 


0 

0 

— ljv 

0 / 


the explicit form of such a matrix is 

b c d \ 

a T -d T -ct 

-d T —a T 6 t 

—c b —a ) 


H = 


t a . 


c 

\d 


1 _ 


where all entries are complex N x N matrices and a, d (b, c ) are hermitian (skew). 

The Gaussian random-matrix ensemble on P is defined by the Gaussian measure dp{H) with vanishing first moment, 
and second moment 


Tr(AH)Tv{BH)dp{H) 


AN 


Tr (AB - ACE T C _1 + ATB t T 


A(CT)B(CT)~ 1 ) . 


Given the auxiliary space W := C 1 ! 1 ® C 2 ® C 2 <S> C™, we impose on ip £ Hom>(W / , V), ip £ Horn? (T, W) the linear 
conditions 


ip = 'yip T C 1 , ip = —Cip T 7 1 , 
ip = rip T T~ 1 , xp = Tip 1 t^ 1 , 

with some invertible orthogonal elements 7 , r of Endo(TT). Consistency requires 7 2 = —a = r 2 and 77 - + T 7 = 0. A 
possible choice is 


7 = ( Ebb ® icr y I 2 + Eff ® o x <E> ox) ® 1„, 
r = (Ebb ® cr z <8> zcr y + E F f <8> I 2 ® cr x ) ® l n - 

Because this differs from class Cl only by the exchange of the bosonic and fermionic subspaces, the following devel¬ 
opment closely parallels that for Cl, and we arrive at another variant of Theorem 3.3. 

The large-TV limit is dominated by a pair of complex-analytic saddle-point supermanifolds, each being isomorphic 
to Osp(2n|2?r). (The reason why there are two is that 0(2?r,C) has two connected components.) The first one is the 
orbit under Ad(GA) of Q 0 = *u1b|f ® cr z <8> I 2 ® 1 n , and the second one is the orbit of 

Qi = iv I Ebb 0 <r z 0 h 0 l n + -Eff 0 (I 2 0 o x 0 Eu + a z 0 1 2 0 ^ Eu) 

V i =2 

Both saddle-point supermanifolds are Riemannian symmetric superspaces of type C|E with dimensionality m = 2 n 
(Table 2). 



D. Class AIII 


The tangent space at the origin of U (p, q)/U(p) x U(g) consists of the matrices of the form 


H = 


0 Z 
0 


where Z is complex and has dimension p x q. Such matrices are equivalently described by H' = H = —VHP -1 where 
V = diag(l p , —1 q ). For simplicity, we will consider only the case p = q (the general case has not yet been analyzed in 
the present formalism). The Gaussian ensemble of random matrices H is taken to have second moment 

r v 2 

/ Tr(AH)Tr(BH)dp(H) = — Tr (AB - AVBV~ X ) . 

The physical space is V = C 2 <S> C p , and the auxiliary space is IT = C 1 ! 1 <g> C 2 <g> C n . The definition of u> is unchanged 
from class C. To implement the symmetry condition ipip = —VipipV^ 1 we set 


21 



ip = iVipw , ip = iiripV 1 

where ir = 1 B | F 0 ia v 0 1„. This choice is consistent with the relation ipB = /3ipB w hich ensures convergence of the 
{ip, ^(-integration. The auxiliary variable Q ranges over the complex-analytic superspace 

Q = {Q € EndA(W)|<5 = -ttQtt" 1 }, 

and the normalizer of Q in Gl(IT) is: 

G a = {g g Gl(W)|ff = Trgn- 1 } ~ Gl(n|n) x Gl(n|n). 

For the integration domain IA in the FF-sector we again take the antihermitian matrices in Q FF . In the BB-sector we 
set 


M = {Y G End c (W B )|F = ttFtt- 1 = -0Y0- 1 = yt}, 

V ± = {X G Endc(Ws)|X = -ttXtt- 1 = ±/3X(3 ~ 1 = T^}. 

The treatment of Sec. EE then goes through as before, leading again to Theorem 3.3. 

There is a single dominant saddle-point supermanifold, which is the Ad(GA)-orbit of Qo = u>l B i F <S> cr z 0 l n and 
is diffeomorphic to G\/Ha ~ Gl(n|n). The integration domain Mb x M f is given by M B ~ Gl(n,C)/U(n) and 
M F = U(n). The invariant Berezin measure of this Riemannian symmetric superspace of type A\A was discussed for 
n = 1 in Example 2.4. 


E. Class BDI 

The form of the random-matrix Hamiltonian H for class BDI can be obtained from the preceding case by taking 
the p x q matrix Z to be real. Put in formulas, H is subject to W = H = H T = We again make the 

restriction to p = q. The basic correlation law of the Gaussian ensemble is 

r v 2 

/ Tr{AH)Tr{BH)dg{H) = — Tr (A{B + B T ) - AV{B + B T )V~ 1 ) . 

To accommodate the extra symmetry H = H T , auxiliary space is extended to W = C 1 ! 1 0 C 2 <g> C 2 ® C". The 
symmetry conditions on ip, ip are 

ip = iVipn -1 , ip = mipV ~ 1 ; ip = ip T T~ 1 , ip = rip T , 

where ir = l B i F (g> ia y 0 1-2 0 In and r = {Ebb 0 I 2 0 o x + Eyy 0 12® i<J y ) 0 I n - The auxiliary integration space 

S = {Qe EndA(W)|Q = —7rQ7r _1 = +t(2 t t“ 1 } 
has symmetry group (or normalizer) 

G a = {g e Gl(W)|g = = rg^r^ 1 } ~ Gl(2n|2n). 

For the integration domain U in the FF-sector we once again take the antihermitian matrices in Q FF . In the BB-sector 
we set 


M = {Y e End c (W B )|y = ttHtt" 1 = = -pYfi- 1 = y f }, 

V ± = {X e Endc(W B )|X = -ttAtt" 1 = +tX t t~ 1 = ±0X0~ 1 = =FX f }. 

The treatment of Sec. EH then goes through wit modifications as in Sec. Evb| . 

There is a single dominant saddle-point supermanifold, which is the Ad(GAj-orbit of Qo = *ul B | F 0 cr 2 0 I 2 0 1 n 
and is diffeomorphic to Ga/Ha — Gl(2n|2n)/Osp(2n|2n). The integration domain Mb x My is given by M B ~ 
Gl(2n, R)/0(2n) and My = U(2n)/Sp(n). This is a Riemannian symmetric superspace of type AI|AII with m = 2 n 
(Table 2). 
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F. Class CII 


The tangent space at the origin of Sp {N, N)/Sp{N) x Sp(iV) (a noncompact symmetric space of type CII) can be 
described by the equations 

= H = - VHV~ l = -TH t T~\ 

where V = a z (g> I 2 g ljv and T = I 2 g icr v g ljv (the physical space is V = C 2 g C 2 g C w ). The explicit form of the 
matrices is 


/° 

0 

a 

b \ 






0 

0 

-b 

a j 






a ^ 

-b T 

0 

0 ’ 






\b* 

a T 

0 

0/ 






( 1 N 

0 

0 

0 \ 

/ 

0 

1 ,V 

0 

0 \ 

0 

1 ,V 

0 

1 and E = 


—ljv 

0 

0 

0 

0 

0 

— In 

o 


0 

0 

0 

ljv 

V 0 

0 

0 

— ljv / 

V 

0 

0 

— ljv 

0 J 


where a and b are complex and have dimension N x TV. The correlation law of the Gaussian random-matrix ensemble 
of type CII is 


Tr(AH)Tr(BH)dn{H) = —Tr(A - VAV~ 1 ){B - TB t T~ 1 ). 


As before, W = C 1 l 1 gC 2 gC 2 g C". The symmetry conditions on ip, ip are 

ip = iVipir -1 , ip = i'iripV~ 1 ; ip = T V’ T r _1 , ip = Tip T T~ l , 


where 7r = 1 B | F g) ia y g 1 2 <g> 1„ and r = (Cbb g I 2 g *cr y + Ep p g 1 2 g a x ) g 1„. This differs from class BDl only by 
the exchange of the bosonic and fermionic subspaces. Once more we arrive at another version of Theorem 3.3. 

There is a only one complex-analytic supermanifold of saddle-points that dominates for N —> 00 . It is isomorphic 
to that for class BDl. The integration domain Mb x M f changes to Mb — U*(2n)/Sp(n) and M F ~ U(2n)/0(2n) 
(not U(2n)/SO(2n)). This is a Riemannian symmetric superspace of type AII|AI with in = 2 n (Table 2). The group 
U*(2n) is defined as the noncompact real subgroup of Gl(2n, C) fixed by g = CgC -1 where C = ia y g l n . 


G. Class A 


This class for n = 1 was used to illustrate our general strategy in Sec. Ill A. Let us now do the case of arbitrary n, 

(^1 ? ■ , /?! , ..., /? n ) — 


/ II Det 

hn{N) fJi 


H - fa 


H - a. 


dn(H). 


The classes treated so far (C, D , Cl, Dill, AIII, BDl, CII) all share one feature, namely the existence of a particle-hole 
type of symmetry (H = —VHV~ 1 or H = —CH T C~ 1 ), which allows to restrict all a* to one half of the complex 
plane. Such a symmetry is absent for the Wigner-Dyson symmetry classes A, Al, and All, which results in a somewhat 
different scenario as it now matters how many at lie above or below the real axis. For definiteness let 

Im Q.i <0 (i = 1,..., ua), Imay >0 (j = + 1,..., n), 

and set hr = n — ua- 

Auxiliary space is taken to be IT = C 1 ! 1 g C". The definition of u> changes to 

n n 

U = I?BB ® ^2, a iEii + Dpp g ^2 PjEjj- 

*= 1 3 =1 


Recall that the imaginary parts of the cti steer the convergence of the {ip, ^-integration. Since ui couples to ip, ip by 
exp —iSTriyw^'i/), convergence forces us to take ipB = ftipB where 
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riA n 

@ y Eji — Ejj. 

i—l j=n A +1 

There are no further constraints on ip, ip, or Q. Thus the complex-analytic auxiliary integration space is Q = EndA(IF), 
and Ga = Gl(IF) ~ Gl(n|n). 

The integration domain for Q in the FF-sector is taken to be the antihermitian matrices U = u(n). In the BB-sector 
we introduce 


g = {x e g i(n, c)\x = -px^p- 1 }, ic = {x £ g\x = pxp- 1 }. 

The Lie algebra Q is a noncompact real form u(n^,n/;) of gl(n, C), and K, = u(tia) © u (n R ) is a maximal compact 
subalgebra. The space M is defined by the Cartan decomposition Q = /C © M. The integration domain for Qbb is 
taken to be i(pb{fc x Ad) where (pb{X, Y) = b(X + e ad ( Y )/3) (b > 0). This gives Theorem 3.3. 

By simple power counting, the limit N —> oo is again dominated by a single complex-analytic saddle-point super¬ 
manifold, which is the Ad(GA)-orbit of Qo = w1b|f©/?- The stability group Ha of Qo is Ha = G\(nA\riA) x Gl(n_R|nij), 
so 


Ad(G A )Qo ^ Ga/Ha = Gl(n|n)/Gl(ruM x Gl(n*|n*). 

The intersection of Ad(GA)Qo with i<Pv{fc x Ad) x U is M B x M F where M B ~ \J(ua, n R )/V(nA) x U (n R ) and 
M f ~ \J (ua + Hu)/\J( ua) x U(n R ). This is a Riemannian symmetric superspace of type AIII|AIII with m i = n i = Ua 
and m 2 = n -2 = n R (see Table 2). 


H. Class AI 

The tangent space of XJ(N)/0(N) is the same as (i times) the real symmetric matrices H' — H — H T . It differs 
from the tangent space of SU(AT)/SO(iV), a symmetric space of type AI, in an inessential way (just remove the 
multiples of the unit matrix). The Gaussian ensemble over the real symmetric matrices has its second moment given 

by 


r v 2 

/ Ti(AH)Tr{BH)dg(H) = —Tt(AB + AB t ). 

This ensemble is related to type A in the same way that type GI is related to type G. 

To implement the symmetry H = H T we set W = C 1 ! 1 © C 2 © C" and require ip = ip T r~ 1 , ip = rip T where 
t = (-Ebb © <J X + -E’ff © i<?y) © 1 n - The auxiliary integration space 

Q = {<5 G EndA(kF)|<2 = t(3 t t -1 } 

has the symmetry group 

G a = {g £ Gl{W)\g = rg^ lT r _1 } ~ Osp(2n|2?i). 

The intersection U of the FF-sector Qff with the antihermitian matrices is given by sp(n) ®U = u(2n). In the 
BB-sector we put 


M = {Y £ End c (W B )|F = - tY t t~ 1 = - PYP= Y t}, 

V ± = {X £ Endc(IT B )|A = +tX t t ~ 1 = ±pxp~ l = ®X^}, 

which leads to yet another version of Theorem 3.3. 

The large -N limit is controlled by a single complex-analytic saddle-point supernranifold Ad(GA)<3o — Ga/Ha where 
H a ^ OsTp{2n A \2n A ) x Osp(2n R \2n R ) is the stability group of Q 0 = w1 B |f © e u - J2j= nA +i E ii) • Tlie inter¬ 

section of Ad(GA)Qo with the integration domain (p v (V + xM)xU is Mb x Mf where Mb — SO(2n A , 2n R )/SO(2n A ) x 
SO(2 n R ) and M f ~ Sp(riA + n R )/Sp(nA) x Sp(tir). This is a Riemannian symmetric superspace of type BDl\Cll 
(Table 2) with mi = 2ni = 2n A and m 2 = 2ri2 = 2 n R . 
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I. Class All 


Finally, the tangent space of U(21V)/Sp(7V) (a symmetric space of type All, except for the substitution SU(27V) —> 
U(27V)) can be described as (i times) the subspace of Encl(C 2 ®C JV ) fixed by the linear equations H 1 = H = TH t T ~~ 1 , 
T = i<jy (g) lj\r. The explicit matrix form of H is 


H = 



b 


where b is skew and a is hermitian. 

The conditions if) = Ti/) T r _1 and = rip^-T -1 are mutually consistent if, say, r = (Ebb ® icr y + Epf ® &x) l n - 
The rest of the manipulations leading up to Theorem 3.3 are the same as for class AI, except for the exchange of the 
bosonic and fermionic subspaces (tb <-> tf). The large- N limit is controlled by a single saddle-point supermanifold 
(Ga/H\,Mb x Afp) where 


Ga/H\ = Osp(2n|2n)/Osp(2nA|2nA) x Ospi2n R \2n R ), 

Af B = Sp(nA,n R )/Sp(n A ) x Sp(n_R), 

Afp = SO(2rc,4 + 2n R )/SO(2riA) x SO(2n^), 

which is a Riemannian symmetric superspace of type CII\BDl (Table 2) with mi = 2ni = 2 n A and m 2 = 2?r 2 = 2,n R . 


V. SUMMARY 

When Dyson realized |ff2| that the random-matrix ensembles he had introduced were based on the symmetric 
spaces of type A, AI and All, he wrote: “The proof of [the] Theorem ... is a mere verification. It would be highly 
desirable to find a more illuminating proof, in which the appearance of the [final result] might be related directly to the 
structure of the symmetric space...”. The advent of the supersymmetry method of Efetov and others has improved 
the situation lamented by Dyson. The present work takes the Gaussian random-matrix ensembles defined over 
Cartan’s large families of symmetric spaces and, going to the limit of large matrix dimension, expresses their spectral 
correlation functions as integrals over the corresponding Riemannian symmetric superspaces. These correspondences 
are summarized in Table 3. 


RMT 

comments 

RSS 

dimensions 

A 

Wigner-Dyson (GUE) 

AIII AIII 

mi = ni = n A , m 2 = n 2 = n R 

AI 

Wigner-Dyson (GOE) 

BDl\Gll 

= 2ni = 2n A , m 2 = 2n 2 = 2n/j 

All 

Wigner-Dyson (GSE) 

Cll\BDl 

mi = 2ni = 2n A , m 2 = 2n 2 = 2 n R 

AIII (p = q) 

chiral GUE 

A A 

rn = n 

BDl [p = q) 

chiral GOE 

AI All 

m = n 

cil ip = q) 

chiral GSE 

All AI 

m = n 

C 

NS 

Dill Cl 

m = n 

Cl 

NS 

DC 

m = 2n 

D 

NS 

Cl | Dill 

m = n 

Dill 

NS 

CD 

m = 2 n 


Table 3: The symmetric-space based random-matrix theories of the first column map onto the Riemannian symmetric 
superspaces listed in the third column. The notation for the dimensions is taken from Table 2. 

The Riemannian symmetric superspaces that appear in Table 3 all have superdimension ( [p , q) with p = q. We say 
that they are “perfectly graded” or “supersymmetric”. An interesting question for future mathematical research is 
whether our procedure can be optimized by reducing it to a computation involving no more than the root system of 
the symmetric space, thereby obviating the space- and time-consuming need to distinguish cases. (Although I have 
treated all ten cases separately, it is possible, following Efetov |j), to shorten the derivation by starting from a large 
“master ensemble” of highest symmetry and then reducing it by the addition of symmetry-breaking terms. I chose 
not to follow this route as it involves handling large tensor products, which makes the computations less transparent 
and the identification of the spaces involved more difficult.) 


25 










The great strength of the supersymmetry method, as compared to other methods of mesoscopic physics, stems from 
the fact that it easily extends beyond the universal random-matrix limit to diffusive and localized systems. What one 
obtains for these more general systems are field theories of the nonlinear er model type, with fields that take values in a 
Riemannian symmetric superspace. The method also extends beyond spectral correlations and allows the calculation 
of wave function statistics and of transport coefficients such as the electrical conductance (see the literature cited in 
the introduction). 

Let me end on a provocative note. Mathematicians and mathematical physicists working on supermanifold the¬ 
ory have taken much guidance from developments in such esoteric subjects as supergravity and superstring theory. 
Wouldn’t it be just as worth while to investigate the beautiful structures outlined in the present paper, whose physical 
basis is not speculative but firmly established, and which are of direct relevance to experiments that are currently 
being performed in physics laboratories all over the world? 

This research was supported in part by the National Science Foundation under Grant No. PHY94-07194. 
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